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We consider quantum many body systems as quantum channels and study the relation between 
the transfer quality and the size of the spectral gap between the system's ground and excited states. 
In our setup two ancillas are weakly coupled to the quantum many body system at different sites, and 
we study the propagation of an excitation and quantum information from one ancilla to the other. 
We observe two different scenarios: a slow, but perfect transfer if the gap large and a fast, but un- 
complete transfer otherwise. We provide a numerical and analytical approach as well as a simplified 
physical model explaining our findings. Our results relate the potential of spin chains acting as 
quantum channels to the concept of quantum phase transitions and offer a different approach to the 
' characterisation of these. 
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I. INTRODUCTION 

> : 

ly-j i Most quantum information processing tasks require at some stage the transfer of quantum states between two 
00 ! quantum systems such as atoms or ions which are located a different positions in space. For transfer over long 
distances, photons sent through optical fibres seem promising. However the interactions between photons and the 
stationary systems, e.g. atoms, are weak and need to be controlled with high precision for transferring the state onto 
the photon and vice versa. Finding alternative methods and carriers is thus of considerable interest, in particular 
for transfer over short distances. Here, using condensed matter systems, e.g. a piece of solid, seems very appealing. 
Therefore the possibilities of transferring quantum information with strongly coupled quantum many body systems 
such as spin chains have been studied in some detail in recent years and several scenarios showing close to perfect 
^p-t state transfer have been found pdllT^|. 

One question arises naturally in this context: How do the properties of the employed many body system relate to 
the transfer quality and speed? A key property in this context is whether those systems feature an energetic gap 
between their ground state and excited states. Most interestingly, this gap vanishes at the critical points of quantum 
^j*. phase transitions where, at zero temperature, the ground state and an excited state exchange their roles as a 
!LJ ' parameter in the Hamiltonian, such as a magnetic held, is varied. 
. £h [ As with classical phase transitions, quantum phase transitions are usually analysed in terms of the scaling behaviour 
of equilibrium properties, where a diverging correlation length is indicative of a critical point Q, 0, Q ■ An analogous 
scaling phenomenon was recently also found for the entanglement properties of a spin chain in the vicinity of a 
quantum phase transition 0, y| Q . Motivated by these findings and the recent experimental observation of the Mott 
quantum phase transition in the well-controlled environment of an optical lattice 0, the dynamical entanglement 
properties of quantum many body systems undergoing a quantum phase transition are receiving increasing attention. 
For example, one recent approach Q was concerned with the dynamics of bipartite entanglement in spin chains 
resulting from an initial perturbation, while another studied the entanglement of two spins that are globally coupled 
to a quantum critical system 0- On another level, the Zurek-Kibble mechanism for classical phase transitions was 
recently generalised to its quantum analogue, further deepening our insight into the dynamics of quantum phase 
transitions [Toj . 

In this article we study the relation between the size of the spectral gap of a quantum many body system and its 
capacity to transfer quantum information. Specifically, we study the transfer of quantum states for two examples of 
linear chains of interacting quantum systems. We employ newly developed matrix product state techniques |l3| to 
simulate numerically the dynamics of spin chains exhibiting a quantum phase transition. Then we proceed to study 
a harmonic chain where we may choose the on-site potential such that the energy gap above the unique ground state 
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FIG. 1: The topology for the spin model considered in the numerical simulations. S labels the sender and R the receiver 
ancilla, while ms and ma label the spins of the chain where S and R couple to. 



vanishes. The latter of the two models allows us to obtain a better understanding of the relevant physics since it 
permits an analytical study in terms of master equations and the verification of the validity of the assumptions inherent 
in the master equation by numerically simulating the dynamics of the harmonic chain with up to 1400 constituents 

m 

We find that the transfer properties crucially depend on the energy gap between the ground state and the lowest 
excited states, but does not significantly dependent on the detailed structure of the Hamiltonian. In particular the 
transfer The characteristics of the state transfer through such systems may therefore be used to detect the critical 
point experimentally. 



II. SPIN CHAINS 



We begin by considering a 1-D chain of spins with nearest neighbour interactions and open boundary conditions. 
The Hamiltonian of our model reads 

N N-l 

ffchain = B Of + (■'•«■ I + J V°i + J *°i°i+l) , C 1 ) 

i=l i=l 

where N is the number of spins, B is an applied magnetic field and J x , J y and J z the interaction between neighboring 
spins. Furthermore, two ancillas (named S for "sender" and R for "receiver" ) couple to the chain at spins ms and 
rriR, which are near the centre of the chain in order to avoid boundary effects . The complete Hamiltonian is thus 
given by 

H = i/ c hain + B a (a z s + a z R ) + J a (of< s + o£< J . (2) 

B a > is the Zeeman splitting of the ancillas, which might differ from B, and J a > is the coupling of the ancillas 
to the chain, which is taken to be weak, i.e. J a <C (B, J X) J y , J z ). Figure ^ shows the topology of the model. 

Initially, the chain is assumed to be in the ground state, |0 c hain)j of the Hamiltonian while the sender is spin 
up and the receiver is spin down. Hence, the initial state of the total system is 

|*(0)> = | Ts,|fl,0 ch ain) ■ (3) 

We simulate the dynamics of our system numerically, making use of the recently introduced matrix product states 
|l3j. We use matrices of dimension 10 x 10. To test the accuracy of our simulations, we verified whether the results 
where stable with respect to variations of the matrix dimension and the size of the timesteps. Furthermore, we tested 
whether the energy of the total system was conserved. Since the matrix product approximation can only be efficient 
if the considered system obeys a "entropy and area law" ^4|, which is not necessarily true at quantum critical points, 
our simulations consider only parameters near, but not exactly on the critical point. 

Figure |21 shows the probability P(tslfl) that "sender" S is in its excited state | ts) and the "receiver" R in its 
ground state | Ir), together with PQsTfl) anci P(IsIr) for a model with N = 100, ms = 45, m^ = 55, B = 1, 
J x = 0.3, J y = J z = 0, B a = 0.64 and J a — 0.05. P{]s]r) is always less than 10~ 4 . The plots show that the 
excitation that was initially located in S oscillates back and forth between S and R. 
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FIG. 2: P(lsln)(t) (blue), P(UU)(t) (green) and P(isU)(t) (red) for B = I, J x = 0.3, J y = J z = 0, B a = 0.64 and 
J a = 0.05 as given by the simulation for the open boundary model with N = 100 spins. S couples to spin 45 and R to spin 55. 
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FIG. 3: P(lsU)(t) (blue), P(Ts!fl)(t) (green) and P(lsU)(t) (red) for B = 1, J x = 0.3, J y = J z = 0, B a = 0.8 and J a = 0.05 
as given by the simulation for the open boundary model with N = 600 spins. S couples to spin 295 and R to spin 305. 



Figure shows P(Tsifl), P(isU) and P(lsln) for a model with N = 600, m s = 295, m R = 305, B = 1, J x = 0.3, 
J y = J 2 = 0, B a = 0.8 and J a — 0.05. Again, P(tsTij) is always less than 10 -4 . For these parameters, the excitation 
is not fully transferred to R, contrary to figure [5] Both, S and R relax to their ground states with the excitation 
only being partially and temporarily transferred to R, even for close-lying spins. Note that the parameters chosen in 
figures |21 and |3 are the same except for B a which in figure [3] is significantly larger than in figure 

The two observed scenarios are rather generic. To demonstrate this, we have done the same simulations for different 
parameters, i.e. for a XXZ-model. The results, shown in figures 01 and clearly agree with our findings for the previous 
coupling parameters. Again B a in figurc|S]is significantly larger than in figure^ while all other parameters are equal. 
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FIG. 4: P(lsla.)(t) (blue), P(UU)(t) (green) and P(lsU)(t) (red) for B = 1, J x = 0.5, J B = 0.2, J 2 = 0.1, B a = 0.04 and 
J a = 0.05 as given by the simulation for the open boundary model with N = 100 spins. S couples to spin 45 and R to spin 55. 
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FIG. 5: P(|sU)(*) (blue), P(UU)(t) (green) and P(ist»)(*) (blue) for B = 1, J x = 0.3, J„ = 0.2, J, = 0.1, B a = 0.2 and 
J a = 0.05 as given by the simulation for the open boundary model with N = 600 spins. S couples to spin 295 and R to spin 
305. 



III. HEURISTIC PHYSICAL PICTURE 

The dramatic difference between the almost perfect transfer scenarios in figures [21 and 0] and the damped scenario 
in figures |3 and [5] has a simple physical explanation. The dynamics we have simulated is given by the Schrodinger 
equation containing the Hamiltonian As a consequence, all moments of the Hamiltonian are conserved, 

(^(t)\H n \y{t)) = (*(0)|ff n |*(0)) = const for any integer n (4) 
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FIG. 6: Sketch of the energy levels of the system. The dots indicate the occupations of the initial state. For this initial state 
only the energy levels between the two horizontal dashed lines are accessible, resulting in almost perfect transfer for the left 
scenario and damping for the right one. 



The initial state 1^(0)) is not an eigenstate of H as given by (0), hence 

(*(o)|iT|*(o))=£|(i?|*(o)>| 2 £, (5) 

E 

where E and \E) are the eigenvalues and eigenstates of H . However since a probability distribution is entirely 
determined by all moments, {^(t)\H n \^(t)) = const for all n implies \{E\^{t))\ 2 = |(£;|*(0))| 2 for all \E). In other 
words the whole probability distribution given by the | (E\^(0)) | 2 is conserved. In our case, it's variance is 

y/(*(t)\H 2 Mt))-m)\HMt)) 2 = J°- (6) 

For the dynamics this means that only those states with an energy expectation value E in the range ('l'(O)l-ff ^(O)) — 
2J a < E < (^(0)1^1^(0)) + 2J a are accessible. Figure El sketches the energy levels of the system we consider. S 
and R are depicted as two level systems, while for the chain there is a unique ground state and a quasi continuous 
band of excited states sketched as the gray area. The dots indicate the initial occupations. The energy range which 
is accessible for the considered initial state lies between the two horizontal dashed lines. If the spectral gap is larger 
than the Zeeman splitting of the ancillas (left plot), there is no accessible excited state of the chain and hence no 
excitations get lost into the chain, which in turn implies the excitation will be almost completely transferred to R. 
If however the spectral gap is smaller than the Zeeman splitting of the ancillas, there are accessible excited states in 
the chain and excitation and hence quantum information get lost. 

In order to obtain a more rigorous justification of this simple picture and to underline the generality of our findings, 
we now turn to a different model for the chain which also features an adjustable energy gap above its unique ground 
state. 



IV. HARMONIC CHAIN 



We consider a harmonic chain with periodic boundary conditions (see fig. [7J described by 



1 N 

tfchain = ^ E (P 2 3 + ^ 2 fe - 9j + l) 2 + (?) 

with the pj being the momenta and the qj the positions (qjv+i — qi). In this case, the two ancillas are harmonic 
oscillators that couple to oscillators ms and ttir of the chain. The complete Hamiltonian now reads 



Htot = H chain + H ancillas + Hi (8) 
ncillas = g (PS + + P% + ^Ir) ( 9 ) 

Hi = J a (qsq ). (10) 
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FIG. 7: The topology for the harmonic chain model considered in the analytical approach. S labels the sender and R the 
receiver ancilla, while ms and ran label the oscillators of the chain where S and R couple to. 



Since we are only interested in the time evolution of the ancillas, we derive a master equation for the dynamics of 
their reduced density matrix p(t). For weak coupling J a <C (f2, fio), its equation of motion is given by 

-£ = -J o ^Tr chain {[ J ff / (0,[i?/(s),|0)(0|®a( S )]]} , (11) 

where a(t) and Hj(t) are the density matrix of the ancillas and the interaction between ancillas and chain in the 
interaction picture, respectively: H = Hq + Hi with Hq = i/ c ham + -^ancillas, Hj(t) = exp(iHQt)Hj exp(— iH^t) and 
<j{t) = cxp(iH t)pexp(—iH t). Tr chain is the trace over the degrees of freedom of the chain and |0) denotes the ground 
state of the chain J7J. The right hand side of eq. Ijllfl is an expansion in the coupling strength J a up to second order, 
which is a good approximation if the integral approaches a constant value for t > t*, where J a t* <C 1. Since a only 
changes significantly on time scales t ~ J^ 1 ^> t* , the approximation a(s) « a(t) can be used. Performing the trace 
on the rhs of (|11J) yields 



f (Xr + (X - X 1 )5 jl ) ({a.ala} - 2{a j( ra\))) , (12) 



dt 

j,l=S,R 



where as and an are the annihilation operators of S and R, respectively: qj = {a j+Oj) / ■s/2u} andpj = —i(aj—Oj)-\/u)/2 
for j = S, R. [•, •] and {■, ■} denote commutators and anti-commutators. On the rhs of the above equation, we neglected 
terms which contain two annihilation or two creation operators since they oscillate at high frequencies. The validity 
of this approximation can later be confirmed from the exact numerical solution. The coefficients read 

X Q = Re(C+ sms + C" sms )/2 W , (13) 

X x = Re(C+ smR + C- smR )/2oj , (14) 

Y Q = Im(C+ sms + C msms )/2uj and (15) 

Fi = lm(C+ smR + C- sm J/2a; (16) 

with given by 

CtM= f'dsiOMtMsMe^-^, (17) 
Jo 

where k, I = ms, rnp. Due to the symmetries of the model, the only depend on \k — l\, implying C^ sms = C^ RTnR 
and Cm smR = ^m R m s - Eq. I|12|) is a good approximation whenever 

« C% = const, for t -C J~ x . (18) 

Since the (t) do not depend on J a themselves, there is always a sufficiently small J a such that (|18(1 holds, provided 
lim^oo CuAt) exists. 
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The harmonic chain can be diagonalised via a Fourier transform 12]. In the limit of an infinitely long chain, 
N — > oo, its dispersion relation is 

k 

u? k = 4ft 2 sin 2 - + 0,1, -7r < k < n , (19) 



and the correlation functions read 



(0\ qj (t) qi (s)\0) = -L r dk^ 1 cos((j -Ofeje-*-*^ . (20) 
27r Jo 

These expressions show, that indeed all lim^oo (i) exist except for the case where u> = fto = 0. As in master 
equations for system bath models, we now insert the asymptotic expressions 

C% = Hm = / (0\q k (t) qi ( S )\0) e^-^ (21) 

into eq. (|12fl . This replacement assumes that all internal dynamics of the chain happens on much shorter time scales 
than the dynamics caused by the interaction of the ancillas with the chain. Furthermore, it does not treat the initial 

evolution for short times with full accuracy since liirii_>o C^i (t) — 0(^ C kt ). The obtained master equation is thus 
valid in a regime where the couplings J a are weak enough such that the time it takes for an excitation to travel from S 
to R is completely determined by J Q , i.e. by the time it takes to be transferred into and from the chain. Consequently, 
the speed of sound of the chain is no longer resolved, and differences in the distance between S and R do not matter. 

From eq. 1|12[) we find the following solution for the expectation values of the occupation numbers of S and R, 
ns = Tr(agascr) and tir = Tr(a R aua): 

n s(t) \ _ ( \ <.)„•) i, ... i ...,.,•> /-',„;,» „v„f_-> n (22) 



nR(t) 



| = (A + cosh(2 Jlx x t) ± A_ cos(2 J 2 ?/it)) exp (-2 J 2 x t) 



Here, A+ = " s(0) +" h(0) , A_ = " s(0) ~" fl(0) , x = lim^oo X ,x x = lim^^ X x and y x = lim t ^oo Yl. Note that > 
and Xq > \xi\. 

Inserting the correlation functions into eq. I|21|l and using the relation 

dx f(x) ^ dr e~ lxr = tt/(0) - iV [ dx i& , (23) 

> JO J-oo x 

where V denotes the principal value of the subsequent integral, one sees that Xq and x\ are only non-zero if lu > lu^ 
for at least one mode k, that is if our initial state is in resonance with (i.e. has the same energy expectation value 
as) states where both ancillas are in their ground states and the chain is in one of its lowest-lying excited states (c.f. 
figure EJ). The dispersion relation shows that this only happens for uj > fto- As in figures [21 and [3] or 0] and [SI we thus 
observe two different scenarios: 

If uj < and therefore xo = x\ = 0, the excitation that is initially in S oscillates back and forth between S and 
R at a frequency 2J 2 j/i, i.e. 



n s (t) 
nit(t) 



A+ ± A- cos(2 Jl Vl t) . (24) 



Note in particular that the excitation is entirely transferred to R at times t n — n(ir/ J 2 yi); n = 1,2, . . . . The solution 
(|24|l is plotted in figure [S] for a harmonic chain and ancillas with f2 = 1, u) = 0.5, J a = 0.05, \ms ~ itir] — 9 and 

n = 0.7. 

Figure [HI shows the frequencies 2J 2 yi of the excitation's oscillations between S and R for cases where ui < fto, 
for fl = 1, J a = 0.05 and u> = 0.35 as a function of Hq. As fto — uj decreases, the transfer becomes faster and the 
oscillation frequency increases. 

If, on the other hand, u> > Qo, the chain acts similarly to a bath. Here, £0 ^ 0, x\ ^ 0, and both ancillas relax 
into their ground state transferring their energy into the chain. During this process, however, a fraction of the energy 
initially located in S appears momentarily in R before it is finally damped into the chain. The maximal excitation 
of the receiver throughout the entire evolution depends on the distance \ms — Wr|- For a given initial energy in 
S, a narrow range of the excitation spectrum is relevant for the dynamics. The relation of the wavelength of these 
excitations to the distance \m,g — mji\ determines the maximal transferred fraction of the excitation. Figure ITUI shows 
the solution [23 for a harmonic chain and ancillas with CI = 1, uj = 0.5, J a = 0.05, \ms — ma\ = 9 and ftg = 0.2. 
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FIG. 10: The solution for a 
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FIG. 11: The reference spin C is initially maximally entangled with the sender S, left plot. If an excitation, which was initially 
in S gets perfectly transferred to R, the entanglement will then be shared between R and C, right plot. 



One might try to derive the same type of master equation for the spin chain (JTJ . However, for finite distances no exact 
expression for the time dependent correlation functions is known [lfj • This is due to the fact that the subspaces of odd 
and of even number of fermions cannot be diagonalised simultaneously. An attempt of an approximation restricted to 
only one subspace led to reasonable results for some parameter values but occasionally produced unphysical solutions 
which grew exponentially in time. Therefore, such a master equation approach cannot be considered reliable for our 
spin chains and was avoided. 

To confirm the validity of the master equation approach for the harmonic chain, we compared it to results of a 
numerical simulation of a chain with 1400 oscillators. Since the Hamiltonians of harmonic oscillators and harmonic 
chains are quadratic in the position and momentum operators, Gaussian states (states with a Gaussian Wigner 
function) remain Gaussian throughout the time evolution. For these states the complete dynamics can thus be 
obtained by only considering the evolution of the covariance matrix (see 0] for details). We found good agreement 
between our analytical and numerical solutions, with the relative errors being less than 5%. 

V. QUANTUM INFORMATION TRANSFER 

The observed effects may also be formulated in quantum information language. In this way, one obtains statements 
on the average fidelity achieved for arbitra ry i nput states (subspace fidelity) or the transfer of entanglement (entan- 
glement fidelity), which are closely related For the present setup, suppose there is an additional control spin C 
which does not couple to the rest of the system, but is initially maximally entangled with S, see figure ITT1 A possible 
initial state is 

l*(0)) = (I Tc, Ts)liJ)0) + I Ids, in., 0»/V2. (25) 

The transfer of the entanglement across the chain may now be analysed by consid- 
ering the entanglement between R and C as a function of time. Since we assume 
J a <C (B, J x , J y , J z ), the state | jcs Is, Ifl: 0) is by virtue of energy conservation approximately stationary, 
while the evolution of | Td ts> -liii 0) is the same as above (modulo a phase). Since states with more than one 
excitation are energetically not accessible, the logarithmic negativity 0] for the reduced density matrix of R and C 
can be expressed approximately as 

E N ^\og 2 {P{isU){t) + l) ■ (26) 
Hence, in this model P(Is1r) ~ 1 implies that the entanglement has been transferred perfectly, too. 

VI. CONCLUSIONS 

In conclusion, we have considered entanglement and excitation transfer through strongly coupled quantum many 
body systems. In particular we have studied the dependence of the transfer quality and speed on the size of the 
spectral gap between the ground and the lowest excited state of the considered system. 
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As a first main result, we find that the quality of transfer, and hence the quantum channel capacity, can be almost 
perfect whenever there is a finite, sufficiently large energy gap above the ground state. This opens up a generic way 
to design good quantum channels by using gapped systems, since the gap ensures high transfer quality irrespective of 
the system's details. 

On approaching quantum critical points, the spectral gap shrinks and the transfer decreases in quality, but accel- 
erates. This second main result suggests a possible experimental determination of the energy gap: if one finds in 
an experiment that the energy is not completely transferred from one ancilla to the other, one can infer that the 
energy gap is smaller than the available energy. The bound can be made tighter by lowering the energy available in 
the "sender" ancilla. This procedure might in particular be helpful for cold atom systems in optical lattices, where 
standard spectroscopy is not applicable. 

A quantitative study of the scaling of the transfer quality and time in between the two detected scenarios and its 
relation to critical exponents of various quantum phase transition universality classes should be a subject of future 
research. In that way an approach that originated in quantum information considerations might open up a new way 
to characterise and experimentally detect quantum phase transitions. 
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